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Abstract 

The charged lepton flavour violating (LFV) radiative decays, fi ^ e + j, t ^ ^ + j 
and T — e + 7 are investigated in a class of supersymmetric A4 models with three heavy 
right-handed (RH) Majorana neutrinos, in which the lepton (neutrino) mixing is predicted 
to leading order (LO) to be tri-bimaximal. The light neutrino masses are generated via 
the type I see-saw mechanism. The analysis is done within the framework of the minimal 
supergravity (mSUGRA) scenario, which provides flavour universal boundary conditions at 
the scale of grand unification AIx ~ 2 x 10^® GeV. Detailed predictions for the rates of the 
three LFV decays are obtained in two explicit realisations of the A4 models due to Altarelli 
and Feruglio and Altarelli and Meloni, respectively. 
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1 Introduction 



Neutrino experiments revealed a totally different mixing pattern in the lepton sector compared 
to the one observed in the quark sector and encoded in the CKM matrix. Indeed, two lepton 
(neutrino) mixing angles are large [1], 

sin2 012 = 0.304;^°:°^^ and sin^ ^23 = 0.50+°;^^ (3a), (1.1) 

while the third one, the CHOOZ angle, is limited by the upper bound [Ij: 

sin^ % < 0.056 {3a) . (1.2) 

The mixing observed in the lepton sector is remarkably close to the tri-bimaximal (TB) one 
proposed in |2]: 

/ ^/2/3 l/VS \ 
Utb = -1/V6 1/V3 1/V2 . (1.3) 
V -1/^6 -1/V2 J 

In the case of TB mixing, the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) neutrino (lepton) 
mixing matrix U is given, in general, by: 



U = Utb 



/ 1 

e'^ I , (1.4) 
^0 



■■ "31 



e 



where 021 and 031 are two Majorana CP violating phases [3]. The three neutrino mixing angles 
in the standard parametrisation of the PMNS matrix (see, e.g. [1]) in this case have fixed values 
determined by: 

sin^ (9^2^ = 1/3 , sin^ = 1/2 and sin^ 6*^3^ = . (1.5) 

Thus, in the case of exact TB neutrino mixing, the Dirac CP violating phase is not present in 
the PMNS matrix and the CP symmetry can be violated only by the Majorana phases. 

It has been shown in [SHS] that TB mixing can arise in a class of models in which the 
alternating (tetrahedral) group A4 serves as flavour symmetry group 0- The light neutrino 
masses are generated in these models either through dimension-5 operators or through the type 
I see-saw mechanism. In the following we concentrate on models which contain three right- 
handed (RH) neutrinos and employ the latter mechanism for neutrino mass generation. The 
light massive neutrinos are Majorana particles. The A4 symmetry is spontaneously broken 
at some high energy scale by the vacuum expectation values (VEVs) of a set of scalar fields, 
called Savons, which transform trivially under the gauge symmetry group of the model. In this 
approach, the TB mixing appears at leading order (LO). The LO results get corrected by higher 
order terms which can play important role in the phenomenology of the models. Such corrections 
are suppressed, in general, by a factor e with respect to the LO results, where e denotes the ratio 
of a generic flavon VEV and the cutoff scale A of the theory. The scale A is assumed to be close 
to the scale of grand unification. The parameter e typically takes values between a few x 10~^ 
and 0.05. 

Apart from predictions for the mixing angles, in A4 models one can obtain also predictions 
for the light and heavy Majorana neutrino mass spectrum [7lEl[TU] and constraints on the 
Majorana CP violating phases. More specifically, both types of light neutrino mass spectrum - 



^TB mixing can also be derived from models with an 54 flavour symmetry 
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with normal and inverted ordering, are possible, but in both cases the lightest neutrino mass is 
constrained to lie in certain intervals. As a consequence, one has specific testable predictions for 
the magnitudes of the sum of neutrino masses and of the effective Majorana mass, measured in 
neutrinoless double beta ((/3/3)oi^-) decay (see, e.g. [ilfTT]). Recently, it has been shown that the 
observed baryon asymmetry of the Universe can be generated in the A4 models under discussion 
through the leptogenesis mechanism [10] (see also [12]; for earlier discussions see |13]). The 
CP (lepton charge) asymmetry, leading to the requisite baryon asymmetry, is proportional to 
the forth power of the A4 flavour symmetry breaking parameter e [10]. The observed relatively 
small value of the baryon asymmetry [2] is thus to large extent due to the smallness of E| e. 
Furthermore, the sign of one of the fundamental parameters of the model, which determines the 
light and the heavy neutrino mass spectrum, can be uniquely fixed by the requirement that the 
generated baryon asymmetry has the correct sign [TOj. 

Other phenomenological predictions of this type of models for, e.g. the rates of charged 
lepton flavour violating (LFV) radiative decays ii — >• £j + 7, i > j, i,j = 1,2,3, where ii = e, 
£2 = fi, £3 = T, the electric dipole moments (EDMs) and magnetic dipole moments (MDMs) 
of the charged leptons, have been studied using effective field theory methods in [15]. In this 
approach [15], a new physics scale M is assumed to exist at (1-^10) TeV. The effective dimension- 
6 operator mediating the charged LFV radiative decays and generating new contributions to 
charged lepton EDMs and MDMs, is suppressed by this scale M, instead of the high energy 
scale A. Thus, the rates of the LFV decays ii — )• £ j + 7 and the EDM of the electron, d^, can 
have values close to and even above the existing experimental upper limits H. Assuming that the 
flavour structure of the indicated dimension-6 operator is also determined by the ^4 symmetry, 
one finds that its form in flavour space is similar to the one of the charged lepton mass matrix. 
In [12] the dependence of the branching ratios B{£i — )• £j + 7), the EDMs and MDMs on e has 
been analyzed in detail. It was found that the contributions of the new physics to the EDMs 
and MDMs arise at LO in e, whereas the LFV transitions are generated only at next-to- leading 
order (NLO). It was shown that B{£i — >■ £j + 7) scales as e^, independently of the type of 
the decaying lepton. Correspondingly, all charged LFV radiative decays are predicted to have 
similar branching ratios. The existing stringent experimental upper bound on B{fj, — )• e + 7) can 
be satisfied if the new physics scale M > 10 TeVQ These results were shown to be independent 
of the generation mechanism of the light neutrino masses. 

In this paper we follow a different approach and compute the branching ratios of charged 
LFV radiative decays in two ^4 models which are based on the Minimal Supersymmetric ex- 
tension of the Standard Model (MSSM) with three RH Majorana neutrinos (for an extensive 
list of articles in which the charged LFV radiative decays were studied see [ISIIIS])- We choose 
as framework the minimal supergravity (mSUGRA) scenario, which provides flavour universal 
boundary conditions at the scale of grand unification Mx ~ 2 x 10^6 GeVl. The SUSY breaking 
and the sparticle masses are completely specified by the flavour universal mass parameters mo, 

^The necessary additional suppression is provided by washout efTects [lOj . 

^As is well known, in the minimal extension of the Standard Model (SM) with massive neutrinos and neutrino 
mixing, the rates and cross sections of the LFV processes are suppressed by the factors [TS] (see also [17p 
2.2 X 10~^ \ Uejm^U^j\'^ / < 5.2 x 10~**, Mw, ruj and Uij, I = e, fi, being the mass, light neutrino masses 
and elements of the PMNS matrix. This renders the LFV processes unobservable. 

*It was indicated in [TS] that the experimental upper limit on the EDM of the electron leads to the constraint 
M > 80 TeV, if the CP violating phases associated with the new physics at the scale M have generic values of 
order one. 

^As has recently been discussed, in the context of global supersymmetry (SUSY) 20 , the presence of soft SUSY 
breaking terms in the fiavon sector can lead to additional flavour non-universal contributions to the sfermion soft 
masses. These vanish in the limit of universal soft SUSY terms in the fiavon potential. 
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mi/2 and Aq. As is well known, this class of models contains also two low energy quantities 
tan /3 and sign(^) as free parameters. Off-diagonal elements in the slepton mass matrices, which 
can lead to relatively large branching ratios of the LFV decays ii — )• ij + 7, are generated 
through renormalization group (RG) effects associated with the three heavy RH Majorana neu- 
trinos [21]. In the present article the branching ratios B{li — t- ij + 7) are calculated using the 
analytic approximations developed in |22H24j . In this approach B{ii — >• ij + 7) depend only 
on the generated off-diagonal elements of the mass matrix of the left-handed sleptons, which 
are functions of the matrix of the neutrino Yukawa couplings (or of the neutrino Dirac mass 
matrix), of the three RH Majorana neutrino masses, and of the SUSY breaking parameters rriQ, 
mi/2 and ^0 (see below). We discuss the case of rather generic NLO corrections to the neutrino 
Dirac mass matrix. The numerical results we present are obtained for two explicit realizations 
of A4 models - those by Altarelh and Feruglio (AF) [TJ and by Altarelli and Meloni (AM) [8]. 
In contrast to the results found in the effective field theory approach in [K], where the setup 
of the AF model was considered, the branching ratios B(ii — )• ij + 7) depend in the case we 
analyse on the type of the light neutrino mass spectrum, and do not necessarily scale with the 
symmetry breaking parameter e. 

From the analytical study of the decay rates we can infer, for instance, that in the case of 
light neutrino mass spectrum with normal ordering (NO), all branching ratios are predicted to 
be independent of the symmetry breaking parameter e, in contrast to the results found in the 
effective field theory approach in [15j. The suppression of the branching ratios B{ii — )• ij + 7) 
in our approach is thus a consequence of the off-diagonal elements of the slepton mass matrices 
being generated through RG effects and not of the smallness of e. A crucial difference between 
the predictions for B{ii — t- ij -I-7) of the AF and the AM models is related to the fact that in the 
AF model there exists an additional source of suppression of B{ii — )• ij + 7). This suppression 
is due to the specific pattern of the NLO corrections (or equivalently due to the specific form 
of the VEVs of the flavons contributing at NLO to the neutrino Dirac mass matrix) in the AF 
model. As a consequence, the predictions for the branching ratios of /li — )■ e -|- 7 and r — >• e -|- 7 
decays in the case of light neutrino mass spectrum with inverted ordering (10) in the AF model 
differ significantly from those obtained in the AM model. 

In the numerical study we perform we obtain predictions for the ii ij+'j decay branching 
ratios in the AF and AM models for one specific point in the mSUGRA parameter space, which 
is compatible with direct bounds on sparticle masses and the requirement of having a viable 
dark matter (DM) candidate and the necessary amount of DM. Results for other points in the 
mSUGRA parameter space can be easily derived by modifying a certain rescaling function and 
adjusting tan /3, present in the expressions for B{ii — )• ij + 7). We discuss also the possibility of 
having B{fi —7-6-1-7) > 10^^^ and B{t — t- ^-1-7) ~ 10~^ in the case of a relatively light sparticle 
mass spectrum, which is easily accessible to LHC [f]. We also show plots for specific values of 
the parameters of an A4 model with generic NLO corrections, which illustrate aspects of the 
predictions for B{fi — )• e -|- 7) that cannot be seen in the scatter plots. 

The paper is organized as follows: in Section 2 we present the aspects of the two A4 models, 
relevant for our study and discuss the generic structure of the NLO corrections to the neutrino 
Dirac mass matrix. In Section 3 we identify the LO and NLO contributions in B{ii — )• ij -I-7) in 
the AF and AM models. This is done within the framework of mSUGRA and for the two possible 
types - with normal or inverted ordering, of the light neutrino mass spectrum. In Section 4 we 
present results of the numerical calculations of the branching ratios B{ii — )• ij + 7). We briefly 

®For recent calculations of the rates of charged LFV decays and fi — e conversion rates within the mSUGRA 
framework with three RH neutrinos and SC'(10)-inspired fermion mass matrices, in which also the DM constraints 
are satisfied, see e.g. [25] . 
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comment on results for // — e conversion and the decays ii — t- 3ij in Section 5. Finally, Section 
6 contains a summary of the results obtained in the present work and conclusions. 



2 A4 Models 

The A4 models [TjlS] we discuss in this paper share several features: the left-handed lepton 
doublets I and the three RH neutrinos ly^ transform as triplets under A4. In contrast, the right- 
handed charged lepton fields e"^, /u'^ and r'^ are singlets under A4 0. The Majorana mass matrix 
niM of the RH neutrinos is generated through the couplings: 

aCKi^^) + 6 Kz^Vs) (2.6) 

where (• • • ) denotes the contraction to an A4 invariant H and 995 ~ 3 and ^ ~ 1 under A4. Here 
and in the following we adopt the convention for the group A4 as given in [71[5]. The vacuum 
alignment of and ips achieved, e.g. in [71|8], is given by: 

{ips) = vseAil,l,lY and {0 = ueA (2.7) 

where vs and u are assumed to be complex numbers having an absolute value of order one. The 
(real and positive) parameter e is associated with the ratio of a typical VEV of a flavon and the 
cutoff scale A of the theory. The generic size of e is around 0.01. At the end of this section we 
will specify the range of e in greater detail. The matrix tum can be parametrized as: 



mM= { -Z 2Z X-Z \ . (2.8) 




It contains two complex parameters X and Z which are conveniently expressed through their 
ratio a = \3Z/X\, their relative phase <j) = arg(Z) — arg(X) and \X\. The parameter \X\ 
determines the absolute mass scale of the RH neutrinos. The matrix niM is diagonalized by 
Utb so that: 

Utb = Utb ^ with n = diag(e-''^i/^ e-*'^2/2^ g-i^3/2) (2.9) 

leads to 

Utb^mUtb = diag(Mi, M2, M3) , (2.10) 

Mi being the physical RH neutrino masses. It has been shown in [101 that one can set ipi = 
without loss of generality. Then the phases ip2 and (^3 coincide at LO in the expansion parameter 
e with the low energy Majorana phases Q21 and 031 as defined in eq. ()1.4p . 
The neutrino Yukawa couplings in the models considered read: 

yuii^'l)hu (2.11) 
so that the neutrino Dirac mass matrix has the simple form: [f| 



mD = yu{ 1 ] Vu (2.12) 




^In the AF model they transform as the three inequivalent one-dimensional representations 1, 1' and 1", 
whereas in the AM model all three right-handed charged lepton fields transform trivially under A4. This, however, 
is irrelevant for our discussion here. 

*There might exist an additional direct mass term, as in the AM model, A/(i/'^!/'^). However, this term leads 
to the same contribution as the term ^{r^'^i/''). 

^We use the convention in which the RH neutrino fields are on the left-hand side of the mass matrix and the 
left-handed fields are on the right-hand side. 
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where Vu denotes the VEV of the MSSM Higgs doublet h^- We can define the matrix of neutrino 
Yukawa couphngs as: 

= ^ . (2.13) 

Vu 

The hght neutrino mass matrix arises from the type I see-saw mechanism: 

TTii/ = —rri^m'^^m£) . (2-14) 
It is diagonalized by Utb- The light neutrino masses m^, i = 1,2,3, are given by: 

m = 4^ . (2.15) 

Mi ^ ' 

At LO, the charged lepton mass matrix mi is diagonal in these models. In the AF model 
the charged lepton masses are generated by the coupling to the flavon (pT with its alignment 
((^t) oc (1)0,0)* (and the coupling to a Froggatt-Nielsen field), whereas in the AM model they 
appear due to the couplings with the flavons ipT and ^' having the alignments {(px) oc (0, 1, 0)* 
and (^') 7^ 0. Note, in particular, that the mass of the r lepton stems from a non-renormalizable 
coupling: 

y,(r^/^T)/id/A • (2.16) 

Since m; is diagonal at this level, the lepton mixing originates only from the neutrino sector and 
is given by eq. ()1.4p . 

This LO result gets corrected by multi-flavon insertions, as well as by shifts in the VEVs of the 
flavons. As a consequence, the matrices tum, mo and mi receive corrections. Correspondingly, 
the lepton masses and mixings receive relative corrections of order e. For our study of LFV 
decays, the form of the corrections of the neutrino Yukawa couplings is of special interest. 
Instead of discussing these for the two specific models, the AF and the AM models, we give 
here a general parametrization of the form of these corrections. We start by writing down the 
CO- variants for the case i/*^ ~ 3 and I ~ 3 under A^: 

{vn) = ulh+vlh + v^h ~ 1 (2.17) 

{uH)' = u^^h + v^2h + vlh ~ 1' (2.18) 

{uH)" = v^h + u^h + ulh ~ 1" (2.19) 

(2i/f/i - - U2I3 

2ulh -v^h-ulh I ~ 35 (2.20) 
2z/^/2 - i^lh - T^lh 

( ^Ih - v^h 

{un)A = i^^h - ulh I ~ 3a (2.21) 
\ -v^h + T^lh 

where 'is{A) is the (anti-)symmetric triplet in the product 3x3. As one can see, the structure 
of mu at LO coincides with the structure coming from the ^4 invariant. 

We shall discuss first the contributions which arise at the NLO level through multi-flavon 
insertions. We assume that such contributions arise at the level of one flavon insertions and are 
thus suppressed by e relative to the LO result. This is true in the two realizations which we 
study below numerically 0. All NLO contributions which are of the same form as the LO result 



^''One could also imagine models in which such contributions are suppressed stronger by higher powers of e. 
However, in the majority of the models, the NLO corrections are suppressed by t only compared to the LO result. 
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can be simply absorbed into the latter. Contributions which cannot be absorbed give rise to 
NLO terms of the form: 



yl,{uH)'i^"K/K + y\„{vH)"i;'K/K + y^s{vH)s(khu/ K + y'^A^'^l) A<t>K / (2.22) 

where -0' and tp" stand for flavons which transform as 1' and 1" under A4, respectively. Here (j) 
denotes a triplet under and, for simplicity, we assume that there is only one such contribution. 
For (-0') = w' eA, {ip") = w" eA and {(j)) = {xi, X2, x^Y e A (with w' , w" and Xj being complex 
numbers whose absolute value is of order one) we find that these induce matrix structures of 
the type: 



6mD = yvw"e\ 10 Vu + yv:w'e\ 1 I (2.23) 





+ys e -X3 2x2 -xi Vu + yA^ 



Apart from this type of contribution we could, in principle, find contributions arising at the 
relative order e due to the perturbation of the VEVs of the flavons at this relative order, when 
NLO corrections are included into the flavon (super-)potential. However, the coupling from 
which the LO term in eq. (j2.1ip originates is generated at the renormalizable level, i.e. without 
involving a flavon. Thus, the most general NLO corrections to the neutrino Dirac mass matrix, 
5mD, are of the form given in eq. (j2.23p . In explicit models the term 5mo has usually a special 
form. On the one hand, the flavons in triplet representations have a certain alignment, such as 
(1,1,1)*, (1,0,0)*, (0,1,0)* or (0,0,1)*. On the other hand, in such models usually there exist 
two different flavour symmetry breaking sectors which are separated by an additional cyclic 
symmetry. In most cases each of these sectors contains one triplet of flavons. Considering NLO 
corrections arising at the level of one flavon insertions, we expect that at most flelds from one of 
the two flavour symmetry breaking sectors can couple at the NLO level to give rise to corrections 
to the neutrino Dirac mass matrix. Thus, there is only one flavon triplet contributing to dniD 
at this level. In the specific framework of the AF model, the NLO terms are given by the triplet 
flavon ifT with {(Pt) = fTeA(l,0,0)* {vt is complex with \vt\ ~ 0{1)), so that we find: 



5mD=ysVT€[ 0-1 \ vu + yAVTe{ 0-1 \ vu. (2.24) 





In contrast, in the AM model we find that the triplet (ps with {fs) = vs e A(l, 1, 1)* gives rise 
to the NLO terms such that: 

2 -1 -1 \ / -1 1 

6mD = ysVse\ -1 2 -1 Ivu + y'AVs^ 1 -1 ] Vu . (2.25) 
-1-12/ \ -1 1 

Similar to the neutrino Dirac mass matrix, the matrices niM and mi also receive corrections 
at the NLO level through multi-flavon insertions and shifts in the flavon VEVs. These corrections 



^^A contribution from the flavon ^ transforming as a trivial singlet under A4 can be absorbed into the LO 
result, as we have already indicated. 
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generate small off-diagonal elements in the charged lepton mass matrix mi. If the corrections 
are of general type, the matrix VeL satisfying: 

V}^m\mi VeL = diag(mg, m^, m^) , (2.26) 

has the form: 

VeL ~ I -zXe 1 zee I (2.27) 




-ZBe -Zee 

where 'z denotes the complex conjugate of z. The parameters Zi are, in general, complex numbers 
and \zi\ ~ 0{1). The Majorana mass matrix mjvf of the RH neutrinos also gets contributions 
from NLO corrections Suim, so that it is no longer exactly diagonalized by Utb, i-e. we have: 

VrU^b i^M + Stum) UtbVr = diag(Mi, Ms, Ms) , (2.28) 
where Vr is defined by: 

(1 WAe WBe \ 
-WAe 1 wee \ . (2.29) 
-WBe -wee 1 / 

The mass eigenvalues Mj are expected to differ from those calculated at LO, Mj, by relative 
corrections of order e. Also here the complex parameters Wi have absolute values \wi\ ~ 0{\). 
We show in both matrices, VeL and Vr, the leading term in the expansion in e for each matrix 
element. In the two models we discuss in more detail one finds [7l|8] that due to the structure 
of the NLO terms, not all parameters Zi and Wi in VeL and Vr, respectively, are arbitrary: in 
the AF model we have za = zb = ze with no constraints on Wi, while in the AM model Zi are 
not related, but wa = and wc = 0. 

Finally, we comment in a more quantitative way on the size of the expansion parameter e, 
which in turn entails constraints on the possible size of tan /3 = {h^) / {h^) = Vu/v^. The upper 
bound on e comes from the requirement that the discussed NLO corrections to the lepton mixing 
angles do not lead to too large deviations from the experimental best fit values. The strongest 
constraint results from the data on the solar neutrino mixing angle and implies e < 0.05. A 
lower bound on e can be obtained by taking into account the fact that the Yukawa coupling of 
the r lepton should not be too large. As mentioned, a rather generic feature of the models of 
interest is that the r lepton mass is generated through a non-renormalizable operator involving 
one flavon. As a consequence, the following relation holds: 

V V \ 

rrir ~ |yr|e(/id) = \yr\e--^cos /S l^^l^;^^^^ (2-30) 
where v ~ 246 GeV. Taking rrir at the Z mass scale, mriMz) ~ 1-74 GeV [26], we find: 

O.Ol^ly^l-^. (2.31) 

tan p 

Reasonable values for \yT-\ are between 1/3 and 3. Using ly^-l =3 and tan/3 = 2 gives: El 

e K 0.007 . (2.32) 

This is the minimal value of e in this type of models. For e ~ 0.05 one finds that \yr\ = 3 
corresponds to the largest allowed value of tan (3 = 15. All smaller values of tan /3 > 2 are 
possible as well. In the numerical analysis we fix e = 0.04. In this case the corresponding 
allowed range of tan/3 is 2 < tan /3 < 12. 



^^As is well known, tan/3 cannot be too small ]27]. We allow here for the rather low value of tan/3 = 2. 
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3 Charged Lepton Flavour Violating Radiative Decays 



3.1 Basic Formulae 

We calculate the branching ratios of the LFV processes ii — )• + 7 (n-^ > rn£.) using the 
following expression 



B{ii Ej+j) ^ B{ii Ij + Vi + Vj) Bo{mo,m 



1/2) 



E 

k 



2 

2 



tan^/3 (3.1) 



where £i = e, £2 = n and £3 = r. In eq. ()3.ip 1^ is the matrix of neutrino Yukawa couplings, 
computed taking into account all NLO effects in the basis in which the charged lepton and RH 
neutrino mass matrices are diagonal and have positive eigenvalues: 

% = V^nU^sYuVeL (3.2) 

where il. is introduced in eq. ()2.9p and Yi, represents the matrix of neutrino Yukawa couplings 
in the basis in which the superpotential is defined (see eq. (|2.13p ). We consider the neutrino 
Dirac mass matrix m/j, including the generic NLO corrections given in eq. ()2.23p . The unitary 
matrices Vr and VeL are given in eqs. (|2.29p and (|2.27p . respectively. 

According to the mSUGRA scenario we consider, at the scale of grand unification Mx ~ 
2 X 10^^ GeV, the slepton mass matrices are diagonal and universal in flavour and the trilinear 
couplings are proportional to the Yukawa couplings: 

(m|)jj = {ml)ij = {ml)ij = SijUil , (3.3) 
{A,)i, = Ao{Y,)i, , (3.4) 

{Ae)ij = Ao{Ye)ij , Aq = aQlTlQ , (3.5) 

where m? and m| are the left-handed and right-handed charged slepton mass matrices, respec- 
tively, while m? is the right-handed sneutrino soft mass term. The gaugino masses are assumed 
to have a common value at the high scale Mx- 

= = M~g = mi/2 . (3.6) 

The scaling function BQ{mQ,mi/2) contains the dependence on the SUSY breaking parameters: 



Bo{mo,mi/2) 



"L (3 + al)ml 



Glm% 



87r2 



(3.7) 



In eq. ()3.7p . Gp is the Fermi constant and aem ~ 1/137 is the fine structure constant. The SUSY 
mass parameter mg in eq. ()3.7p was obtained by performing a fit to the exact KG results |23]. 
The resulting analytic expression in terms of mo and mi 12 has the form [23]: 

m| PS 0.5mgm^/2 ("^0 + 0-6W'?/2)^ ■ (3-8) 

According to [23] , deviations from the exact RG result can be present in the region of relatively 
large (small) mi/2 and small (large) tuq. In Fig. [1] we show the dependence of BQ{mo,mi/2) 
on mi/2 for fixed values of mo and Aq = 0. We notice that the function Bq and, consequently, 
the branching ratio in eq. ()3.ip . can vary up to three to four orders of magnitude, depending 
on which point (mo, mi/2) of the parameter space is considered. For smaller values of mi/2 we 
have larger LFV branching ratios or, equivalently, the larger is the SUSY mass parameter ms, 
the stronger is the suppression of the predicted branching ratio. Values of vlo 7^ lead to larger 
values of i?o("iO) "11/2) (see eq. (|3.7p ) and to an increase of the branching ratios, see eq. (|3.ip . 
In the numerical analysis we set Aq = 0, unless otherwise specified. 
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Figure 1: The dependence of the scahng function 80(1710,1711/2) (see eq. p.7p ) on for 
Ao = and fixed mo- i) mo = 100 GeV (black, dotted hue), ii) mo = 400 GeV (red, dot-dashed 
hne), in) mo = 700 GeV (green, dashed hne) and iv) mo = 1000 GeV (blue, continuous line). 



3.2 Identifying the Leading Order Contributions in B(£i £j + 7) 



We work for convenience with LFV branching ratios given in eq. ()3.ip , normalized to the partial 
branching ratios of the ji 01 t decays into lighter charged lepton and two neutrinos: 



B'i 



B{ 



+ 7) 



B{1, 

e + 7), B{t - 



Vi 



e + 7) ~ 0.18 -B'(^ 



(3.9) 

e + 7) and B(t — 



We have: B{[i ^ e + 7) B\ix 
// + 7) Ri0.17S'(r^/x + 7) 128J. 

It proves useful to analyze separately the contributions in the LFV branching ratios, which 
are associated with each of the three heavy RH Majorana neutrinos. For this purpose we 
rearrange the terms in eq. ()3.ip in the following way: 



B'( 



with 



ij + 7) 



(yjf;) . . log 



m\ 



m.2 
m\ 



+ (i;t).3(i;)3^.1og 



m-3 
m\ 



(1.5 X 10"^ eV) yi sin^ /3 f« 1.5 x 10~^ eV 



Mx 

where we have fixed = \ and used sin^ /3 ~ 1, which is a good approximation given the fact 
that tan/3 > 2. In eq. (|3.10p . m^, k = 1,2,3, are the LO neutrino masses in the models, 
given by eq. ()2.15p and generated via the type I see-saw mechanism. We neglect contributions 
to mfc which arise from NLO corrections in the superpotential since the branching ratio depends 
only logarithmically on the light neutrino masses and typically all such relative corrections are of 
order e « (0.007-^0.05), as discussed in the previous section. We also do not consider RG effects 
in the calculation of the neutrino masses and mixings. Such corrections can be relevant in the 



(3.10) 
(3.11) 
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e + 7 




{y:^)2^{%u 


( w"rv + w'Vv' - ^2{y\ + Vs) - My A + ys))<^ + 0{e^) 
H + 0{e) 
0{e) 


r — e + 7 




{Y^%)u 

(^'J)32(^;)21 
(^'J)33(^;)31 


( w"yl, + w'rv, + X2{rA " ^5) + ^sCy^^ - ^5) ) ^ + ^(e') 

0{e) 


r ^ + 7 




{Yi%h2 

{Y^h2{%)22 

{Y^Z'iiY.U 


y, ( u;"^^', + w'y\„ + 2 X2 + 2 X3 ) e + ©(e^) 
+ 0{e) 
-hi + 0{e) 



Table 1: Combination of elements of the matrix of neutrino Yukawa couplings, 1^, which enter 
into the expression for the branching ratios of the LFV decay £i — )• + 7 (see eq. (j3.10p ) . The 
expression for the relevant 0{e) terms in {YjY,y)ij {i 7^ j) is also given (see text for details). 



case of quasi-degenerate (QD) light neutrino mass spectrum. These corrections are relatively 
small or negligible if the spectrum is hierarchical, or with partial hierarchy [29J. In the 
models under discussion, the lightest neutrino mass in the case of NO (10) mass spectrum is 
constrained (by the data on the neutrino oscillation parameters, see, e.g. [Sllin]) to lie in the 
interval 3.8 x 10~^ eV < mi < 7 x 10~^eV (0.02 eV < m^). Thus, in the case of 10 spectrum 
we present results for 0.02 eV < < 0.10 eV. As we will see, the predictions for the branching 
ratios of the LFV decays ii — )■ +7 depend, in general, on the type of neutrino mass spectrum. 

Neutrino Mass Spectrum with Normal Ordering 

In the case of NO mass spectrum, the three logarithms in eq. ()3.10p are all positive and 
are of the same order (see Fig. [21 left panel). The dominant contribution to the decay ampli- 
tude depends strongly on the combination of Y^ matrix elements in eq. p.lOp . Note that the 
matrix elements of 1^ take all 0{1) values, except for the (31) entry which typically scales as 
the expansion parameter e. This is due to the presence of the TB mixing matrix Utb in the 
expression for the neutrino Yukawa couplings Yy, eq. (|3.2p . 

In Table [T] we give the order of magnitude in e of the coefficients of the three logarithms, 
{YuYy)ij, iYu)i2{Yi^)2j and (^'J)j3(^i/)3j, which appear in the three branching ratios B'{£i — )• 
£j + 7) of interest. The VEVs of the flavon fields are assumed to be real for simplicity. As Table 
[T]shows, the coefficient of the log(m2/m'i) term in each of the three LFV branching ratios under 
discussion is of order one. The same conclusion is valid for the coefficient of the log(m3/mi) 
term in B'[t — t- + 7). In what concerns the coefficients of the term cx log(mi/m*), they 
always originate from NLO corrections in the superpotential. These coefficients correspond to 
the off-diagonal elements of the hermitian matrix yJy^, in which the rotation matrices, Utb, 
and Vr, associated with the diagonalisation of the RH neutrino Majorana mass term, do not 
appear. At order e, they depend only on the parameters of the neutrino Dirac mass matrix (see 
eqs. (I2T2I1 and ([223])), as reported in Table III 
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Normal Ordering Inverted Ordering 
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0.0040 0.0045 0.0050 0.0055 0.0060 0.0065 o.02 0.04 0.06 0.08 0.10 
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Figure 2: The three different contributions in B'{ii — t- ij + 7), eq. ()3.10p . in the case of Hght 
neutrino mass spectrum with normal (inverted) ordering (left (right) panel): i) log(mi/m*) 
vs mi (772,3) (continuous line), ii) log(m2/mi) vs mi (ms) (dashed line) and in) log(m3/mi) 
vs mi (771,3) (dotted line). The results shown correspond to the best fit values [1] |A77z^| = 
2.40 X 10-3 eV^ and r = Am|/|A777|| = 0.032. 



Taking into account the magnitude of the different terms shown in Table [U we expect, in 
general, that in the case of NO neutrino mass spectrum: 

B'{n ^ e + 7) pa B'{t ^ e + 7) f« Bo{mo, mi/2) 

B'{t fi + j) ^ Bo{mo, 7711/2) 

Prom the analytical estimates, eqs. (|3.12p and (j3.13p , we conclude that in the case of NO mass 
spectrum, B'{t — )■ ^ + 7) is approximately by one order of magnitude larger than B'{t — )• e + 7) 
and B'^jj, — >■ e + 7). 

Neutrino Mass Spectrum with Inverted Ordering 

As can be seen in Fig. [2] (right panel), the term proportional to log (7712/7771) is strongly 
suppressed with respect to the other terms in the case of 10 mass spectrum. This is valid for 
all values of the lightest neutrino mass allowed in the models of interest, 7773 > 0.02 eV. More 
specifically, one has: log (7772/7711) ~ 0.014 (0.003) for 7773 = 0.02 eV (0.1 eV). In the case of 
non-QD neutrino mass spectrum (7773 < 0.1 eV), B'{£i — ij + 7) are determined practically by 
the sum of the terms proportional to log (7711/777*) and log(777,3/r77,i). The second term increases 
as 7n,3 decreases towards the minimal allowed value 7n,3 w 0.02 eV, so that | log(777,3 / 777i)| ~ 1 
(0.1) for 7n,3 = 0.02 eV (0.1 eV). Thus, taking into account the results reported in Tabled! we 
have in the LO approximation in e: 

B'{fi^e + j) « B'{T^e + j) oc 0(6^), (3.14) 
B'{t — )• ^ + 7) oc 



:yl log 



777 2 
7771 



tan^ /3 oc 0.1 l^i, 



(3.12) 



^y,y log 



7712 
7771 



l^Vu log , 

2 \ 7771 



tan^ /3 oc \yi^\ . (3.13) 



-y^ log — 

2 \ 777i 



0.25 \yu\ , for 7713 = 
0.0025 |7/i.|^, forms 



0.02 eV , 
= 0.1 eV 



(3.15) 
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The shown order of magnitude estimates for — t- e + 7) and B'{t — t- e + 7) in eq. ()3.14p 
can be significantly modified by the rather large contribution of the term containing the factor 
log(mi/m^,) ~ (3.5 -^4.5). It follows from Table [1] that for, e.g. e w 0.04, the contribution 
in the LFV branching ratios due to the indicated term can be ~ e log(mi/m^,) ~ 1/5 ~ ^/e 
such that the branching ratios of the decays /x — t- e + 7 and r — )• e + 7 scale as 0{e). For 
771-3 ~ 0.1 eV, B'{fi — )• e + 7) and B'{t —7-6 + 7) can be comparable to the normalized branching 
ratio of r — )• /u + 7 decay, eq. (|3.15p . Indeed, for ms ^ 0.1 eV and e ^ 0.04, owing to the 
interplay between the leading term in the expansion parameter e, log(m3/mi), whose absolute 
value decreases with increasing of m^, and the contribution from log(mi/m^,), B'{t —7-^ + 7) 
scales as few times e. 

Comparing the results for the NO and the 10 neutrino mass spectrum we see that in a model 
with generic NLO corrections to the matrix of neutrino Yukawa couplings, the magnitude of the 
branching ratio B'{t — )• /u+7) practically does not depend on the type of neutrino mass spectrum. 
For e 0.007, B'{fi — )• e + 7) and B'(t — )■ e + 7) in the case of 10 spectrum can be by one order 
of magnitude smaller than in the case of NO spectrum, while if e ~ 0.04, these two branching 
ratios are predicted to be essentially the same for the two types of spectrum. We always have 
(independently of the type of the spectrum and of the value of e) i?'(/i — )• e + 7) i?'(T — )• e + 7). 

In the next section we study numerically the LFV processes in the AF and AM models. 
One important difference between the two models is in the predicted off-diagonal elements of 
the hermitian matrix Y^Yy. In the AM model they are all of ©(e) and originate from the 
NLO corrections to the Dirac mass matrix, eq. (I2.25p . The exact expressions for the matrix 
elements can be derived using Table [T] and setting w' = w" = and (xi, 2:2, 2:3) = (1,1,1)^5. 
In what concerns the AF model, the VEV structure of the flavon fields, w' = w" = and 
(xi,X2,X3) tx (1,0,0), implies that the leading term in the off-diagonal elements of the matrix 
yjf; is of 0(62). rj,j^ig 

receives contributions from the Dirac mass term, see eq. (j2.24p . as well as 
from the charged lepton sector (through Vei, see eq. ()2.27p ). This difference in the e dependence 
of the elements of yJYi, in the two models leads to different predictions for the LFV branching 
ratios for the 10 neutrino mass spectrum. As a consequence, in the AM model the branching 
ratios of the decays /i — )• e + 7 and r — )• e + 7 are up to two orders of magnitude larger than 
those in the AF model. In contrast, for < 0.1 eV, we expect similar results in both models 
for the decay r + 7 since the coefficient of the term proportional to log (777.3/7711) is of order 

We note that in the case of a QD light (heavy) neutrino mass spectrum, 7713 > 0.1 eV, the 
term proportional to log(77T,i/777=^) in eq. (|3.10p gives the dominant contribution and thus the 
magnitude of the non-diagonal elements of Y^Y^, determines the magnitude of the branching 
ratios of the LFV decays. 

The preceding discussion shows that in the models, the LO structure of the matrix of 
neutrino Yukawa couplings 1^, which is determined by Utb, together with the possibility of 
having a heavy RH neutrino mass spectrum with partial hierarchy, leads to LFV decay rates 
scaling as 0{e^). This prediction differs significantly from the one obtained in the effective field 
theory approach. In [[15] the branching ratios of the charged LFV radiative decays were shown 
to scale as in a generic effective field theory framework, and could even be stronger suppressed 
(scaling as e'') in a specific super symmetric scenario. 

Concerning the absolute magnitude of the branching ratios we remark that these are expected 
to be of similar size in both approaches, because the suppression due to (positive) powers of e 
present in the effective field theory approach corresponds in our case to the suppression factor 
associated to the fact that flavour violating soft slepton masses are generated only through RG 
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running. The scales mg and M, which are the relevant scales for charged LFV radiative decays, 
in our approach and in the effective field theory one pTF], respectively, can be related to each 
other. Assuming that the mass scale M arises from one-loop effects of new particles, such as 
SUSY particles, we see that the mass ms of these new particles is identified with M weighted 
with the coupling g of these particles to the charged leptons and divided by the loop factor in. 
Thus, we roughly have: ms ~ g'M/(47r). 

4 Numerical Results 

In this section we report results of the calculations of branching ratios of the LFV decays 
+ T— >-e + 7 and r —)•// + 7. This is done in the form of scatter plots showing the 
correlations between each two of the indicated branching ratios. The calculations are performed 
in the framework of the AF and AM models. The expansion parameter e is set equal to 0.04 in 
the numerical analyses. 

We consider a scenario in which the sparticle mass spectrum is moderately heavy: 

mo = 150GeV, mi/2 = 700 GeV, = GeV, tan/3 = 10. (4.1) 

The parameters in eq. (j4.1|) lead to squark masses between 1.1 TeV and 1.5 TeV, gluino masses 
around 1.6 TeV, and masses of right-handed sleptons are 300 GeV. Thus, these sparticles are ac- 
cessible at LHC. This point in the mSUGRA parameter space belongs to the stau co-annihilation 
region [30fl32j , in which the amount of DM in the Universe can be explained through the lightest 

§ article (LSP). The latter is a bino-like neutralino and has a mass of approximately 280 GeV 
. As has been shown in [33], the stau co-annihilation and the bulk regions are hardly affected, 
if RH neutrinos are included into the mSUGRA context. For the set of parameters in eq. ()4.ip . 
all decay rates scale with the factor Bo^mQ, mi/2) tan^ /3 w 3.8 x 10""^'^. 

The scatter plots are obtained by varying all the 0(1) parameters that enter in the matrix 
of neutrino Yukawa couplings Y^, defined in eq. (|3.2|) . Some of these parameters are equal to 
zero or have a common value. More specifically, in the AM model wa = wc = and in the AF 
model ZA = zb = zc- The NLO corrections to the Dirac mass matrices for the AF and AM 
models are given in eqs. (|2.24p and (|2.25p . respectively. In the calculations of the normalized 
branching ratios B'{£i — )• £j + 7), we set yu = ^ and the absolute values of all the other (complex) 
parameters in 1^ are varied in the interval [0.5, 2]. The corresponding phases are varied between 
and 27r. 

The results obtained for the AF and the AM models and for both the NO and 10 light 
neutrino mass spectrum are presented graphically in Figs. [3] and HI respectively. The scatter 
plots correspond to three values of the lightest neutrino mass: i) mi = 3.8 x 10~^ eV, 5 x 10^^ 
eV and 7 x 10~^ eV (NO spectrum); ii) m^ = 0.02 eV, 0.06 eV and 0.1 eV (10 spectrum). In 
all numerical calculations we neglect the RG effects on neutrino masses and mixings. This is a 
sufficiently good approximation provided the light neutrino mass spectrum is not QD [29]. In 
the class of models we consider, the latter condition is fulfilled for the NO spectrum since 
the lightest neutrino mass mi is constrained to lie in the interval (3.8 7) x 10~^ eV. In the 
case of 10 spectrum, the condition is approximately satisfied for m^ ^0.1 eV. Correspondingly, 
for 10 spectrum we present results for 0.02 eV < m^ ^0.1 eV. 

^^The sparticle masses quoted above have been calculated with ISAJET 7.69 [33] . 
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Figure 3: Correlation between B'{fi — t- e + 7), i?'(T — t- e + 7) and B'{t — t- /i + 7), calculated in the 
AF model. The results shown are obtained for three different values of the lightest neutrino mass 
for both types of neutrino mass spectrum: i) with normal ordering (left panels), mi = 3.8 x 10~^ 
eV (red x), mi = 5 x 10~^ eV (green +) and mi = 7 x 10~^ eV (blue o); ii) with inverted 
ordering (right panels), ms = 0.02 eV (red x ), ms = 0.06 eV (green +) and ms = 0.1 eV (blue o). 
The horizontal dashed line corresponds to the MEGA bound [35], B'[fi —7-6 + 7) — 1-2 x 10^^^. 
The horizontal continuous line corresponds to B'{fi — )• e + 7) = 10^^^, which is the prospective 
sensitivity of the MEG experiment ^36j. 



4.1 Predictions of the AF model 

The results for the AF model are shown in Fig. [3j In the case of NO spectrum (left panels in Fig. 
[3]), the normalized branching ratios i3'(^ — t- e + 7) and i?'(r — >• e + 7), defined in eq. ()3.9p . are 
approximately the same, as the analysis performed in Section 3 suggested. The branching ratios 
are larger for smaller values of the lightest neutrino mass mi , the dominant contribution being 
due to the term oc log(m2/mi) which is a decreasing function of mi (Fig. [21 left panel). The 
same feature is exhibited by the term cx log(m3/mi). The latter is multiplied by a coefficient 
of 0{e). As we have already indicated, the term oc log(mi/m*) in the AF model is suppressed, 
being of 0{e'^), and has a negligible effect on the results. Due to the fact that the coefficient 
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AM: Nominl Ordering 
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Figure 4: The same as in Fig. [3l but for the AM model. 



of the term oc log(m3/mi) in i?'(r —)■ /i + 7) is of order one, the normahzed branching ratio of 
r —>■// + 7 decay is approximately by a factor of ten larger than those of ;U — )• e + 7 and r — )• e + 7 
decays, which is consistent with the analytic estimates given in eqs. (|3.12p and (j3.13p . 

We observe that, for the set of boundary conditions we have chosen, eq. ()4.ip . the MEGA 
upper limit |35j on B{fj. — t- e + 7) is not satisfied for mi = 3.8 x 10~^ eV. This important 
experimental constraint can be satisfied for larger values of the lightest neutrino mass and, in 
particular, for the two other chosen values of mi, mi = 5 x 10^^ eV and mi = 7 X 10"^ eV. 
However, B{fi — )• e + 7) is always larger than lO""*^^ and thus is within the range of sensitivity of 
the MEG experiment [36], B{fi — e + 7) > 10~^^, which is currently taking data. The predicted 
rates of the r LEV radiative decays are always below the current experimental upper bounds ^37j 
as well as below the sensitivity planned to be reached at a SuperB factory [38] . 

In the case of 10 mass spectrum, the predicted — )• e + 7) is always compatible with the 
existing experimental upper limit [35]. In this case the MEG experiment will probe a relatively 
large region of the parameter space of the model. The branching ratios of — >• e + 7 and 
r — )• e + 7 decays are, in general, smaller by up to two orders of magnitude than in the case of a 
neutrino spectrum with NO. As we have explained earlier, this is partly due to the fact that the 
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term oc log(m2/mi), which in the case of NO mass spectrum gives the dominant contribution, 
is strongly suppressed since m2 and mi are nearly equal, m2 = mi, and partly due to the fact 
that the coefficient of the term proportional to log(mi/mH,) is of order e^. This conclusion is 
valid for all allowed values of the lightest neutrino mass, m^ > 0.02 eV (Fig. [2l right panel). In 
contrast to the case of a NO neutrino mass spectrum, the branching ratios of /x — )• e + 7 and 
r — >• e + 7 decays do not show any significant dependence on the lightest neutrino mass, m^. At 
the same time, the r — )• /i + 7 decay branching ratio exhibits a strong dependence on the value 
of 7713. Indeed, it varies by up to two orders of magnitude when m^ is varied from 0.02 eV to 
0.1 eV (Fig. [3l right bottom panel). The magnitude and the behaviour of B'{t —7-^ + 7) as a 
function of 7713 is determined by the term proportional to log(m3/mi) in the right-hand side of 
eq. (|3.1U|) . It has a maximal value for m^ = 0.02 eV and decreases as m^ increases, following 
the decreasing of log(m3/?7ii). As a consequence of the suppression of the coefficient of the 
log(mi/m*) term, the analytic estimates reported in eqs. (j3.14p and (j3.15p are valid. Thus, the 
r — 7- ^ + 7 decay has a branching ratio which, at least for m^ ~ 0.02 eV, is by approximately 
two orders of magnitude larger than those of the two other charged LFV radiative decays. For 
m3 = 0.02 eV we have B'{t — )• /u + 7) ~ lO"^''. Therefore as like in the case of NO spectrum, 
the predicted B'{t — )• /U + 7) for the values of mSUGRA parameters considered is below the 
sensitivity range of the currently planned experiments. 

4.2 Predictions of the AM model 

Our results for the AM model and both types of neutrino mass spectrum are illustrated in Fig. 
m As was discussed above, the main difference with the AF model is in the prediction for the 
coefficient of the term log(mi/m-*). In the AM model this coefficient is of 0(e) for the three 
radiative decays and, therefore, the term oc log(mi/m=K) is not negligible. Obviously, log(mi/m*) 
is a monotonically increasing function of the lightest neutrino mass (Fig. [2]) . Since the coefficient 
of this logarithm is a number with absolute value of order one, for both types of neutrino mass 
spectrum the /i— >-e + 7, r— >-e + 7 and r — t- /i + 7 decay branching ratios exhibit much weaker 
dependence on the lightest neutrino mass compared to the dependence they show in the AF 
model. Most importantly, as a consequence of the contribution due to the term oc log(mi/m*), 
B'{ii — )• e + 7) and B'{t — )• e + 7) in the case of 10 spectrum are predicted to be of the same 
order of magnitude as in the case of NO spectrum (Fig. H]). This is in sharp contrast to the 
predictions of the AF model. 

We find that the predictions for B'{t — )• /^+7) in the cases of NO and 10 spectrum essentially 
do not differ and are similar to those obtained in the AF model. As Fig. H] shows, for both the 
NO and 10 mass spectrum we get B{fi — )■ e + 7) < 1.2 x 10^^^ in roughly half of the parameter 
space explored. At the same time, in practically all the parameter space considered we find that 
B{fi — e + 7) > 10~^^. The r LFV radiative decays are predicted to proceed with rates which 
are below the sensitivity range of the planned experiments. 

4.3 On the Possibihty of Large B{fi — )■ e + 7) > 10~^^ and B{t — )■ + 7) 
10"'^ in the AF Model 

As we have seen, for the point in the mSUGRA parameter space considered the r — )• e + 7 and 
r —)•// + 7 decay branching ratios are predicted to be compatible with the existing experimental 
upper bounds and below the sensitivity of the future planned experiments. However, the decay 
r — )■ ^ + 7 might have a rate within the sensitivity range of the future experiments if the SUSY 
particle masses are smaller (i.e., the effective SUSY mass scale ms, eq. (|3.8p is lower) than those 
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AF: InvertEd Oideriog AF: Inverted Ordering 

iiH3=TO GeV, Mj^j=jm GeV, Ao=70 GeV,tan^lO m„=70 GeV, mi/2=300 GeV, Ao=70 GeV, tan;8=10 
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Figure 5: Left panel: correlation between B'{fi ^ e + 7) and B'{t ^ + in the AF model for 
three different values of the lightest neutrino mass: 772.3 = 0.02 eV (red x), ms = 0.06 eV (green 
+) and 777-3 = 0.1 eV (blue o). The horizontal dashed line shows the current upper bound from 
the MEGA experiment [35], while the continuous line corresponds to the foreseen sensitivity 
of the MEG experiment [36] . The vertical dashed line indicates the possible future bound on 
r —>•// + 7 from a SuperB factory [38j. Right panel: correlation between B'{t — >• ^ + 7) and the 
effective Majorana mass 777ee. The horizontal continuous line shows the prospective reach of a 
SuperB factory. The two dashed vertical lines indicate the expected sensitivity of the GERDA 
II and GERDA III phase [39]. 



resulting from eq. (|4.ip . This possibility can be realized for smaller values (than those we have 
employed) of the mass parameters 7?7o and mi/2- Indeed, consider the following set of values: 

mo = 70 GeV, mi/2 = 300 GeV, Aq = 70 GeV, tan/3 = 10. (4.2) 

For the values given in eq. ()4.2p squarks can be as light as 500 GeV, gluinos have masses of 
approximately 700 GeV and all sleptons have masses smaller than 250 GeV. The LSP providing 
the correct amount of DM in the Universe is bino-like and has a mass of 115 GeV. The parameters 
given in eq. (j4.2p correspond also to a point in the stau co-annihilation region, very close to 
the region excluded by the LEP2 data |30ft32j : the mass of the lightest Higgs boson is near 
114.4 GeV 0. For the indicated values of the SUSY breaking parameters the predicted LEV 
branching ratios are larger than those corresponding to the mSUGRA point in eq. (|4.ip since 
BQ{mQ,mi/2)taLn'^ f3 « 2.3 x 10^*^. As a result, the AM model is strongly disfavored by the 
experimental limit on B{fi — )■ e + 7). In the AF model the latter constraint cannot be satisfied, 
if the neutrino mass spectrum is with NO. In the case of 10 mass spectrum, however, the 
predicted B{fj, — t- e + 7) is compatible with the MEGA bound in nearly half of the region of 
the relevant parameter space and (with the exception of singular specific points) is within the 
sensitivity reach of the MEG experiment. We show in Fig. \5\ left panel, the correlation between 
the normalized branching ratios of the decays — >• e+7 and r — )• /i+7 in the AF model, assuming 
10 light neutrino mass spectrum. The prospective sensitivity of the searches for the r + 7 
decay, which can be reached at a SuperB factory, B{t — )• /i + 7) 10^^ [38j, is also indicated. 

^*A11 masses have been calculated again by using the program IS A JET 7.69 [55] . 
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Generic: Normal Ordering: C = 1.3 Generic: Inverted Ordering: C = 1 



mo=150 GeV, mi/2=700 GeV, Ao=0, tan)S=10 »io=150 GeV, mi/2=700 GeV, Ao=0, tan;8=10 




Figure 6: — t- e + 7) vs rriee for NO (left panel) and 10 (right panel) light neutrino mass 
spectrum calculated for an model with generic NLO corrections, see eq. ()2.23p . Lower and 
upper limits on i?'(^ — )• e + 7) are shown, which can be found by using eq. (|4.3|) for all possible 
combinations of cri,2,3- The horizontal dashed line corresponds to the MEGA bound [35], B'{^ — )• 
e + 7) < 1-2 X 10^^^. The horizontal continuous line corresponds to B'{pi — )• e + 7) = 10"^^, 
which is the prospective sensitivity of the MEG experiment [36j. The results shown correspond 
to the best fit values ^ |Am^| = 2.40 x lO'^ eV^ and r = Am|/|Am^| = 0.032. 



Assuming a scenario in which in the MEG experiment it is found that — )• e + 7) > 10~^^ 
and the SUSY particles with masses, as predicted above, are observed at LHC, we see from 
Fig. O left panel, that B{t — )■ /i + 7) might be detectable at a SuperB factory if the lightest 
neutrino mass w 0.02 eV. For 7713 = 0.02 eV, the (/3/3)o;^-decay effective Majorana mass is 
predicted [10] to lie in the interval rriee ~ (0.018 0.054) eV. Values of rriee in the indicated 
interval might be probed in some of the next generation of (/3/3)oi/-decay experiments (see, 
e.g. |39y40j). In Fig. [Sj right panel, we show the correlation between the normalized branching 
ratio of r — )• ;U + 7 decay and the effective Majorana mass ruee- The relation between niee and 

the lightest light neutrino mass is discussed in [8l|T0l|12]: niee = \J + \ Am\ | 1 2 + e*"^^ | /3, 
where ms and 021 are both functions of one parameter and thus their values are correlated. We 
indicate the prospective sensitivity of the GERDA II and GERDA III phase, mee = 0.09 eV and 
niee = 0.02 eV, respectively [39]. As one can see, with a positive signal of B{t — >• ^ + 7) 10~^ 
at a SuperB factory values of rUee up to niee ~ 0.04 eV can be probed. The sum of neutrino 
masses reads [10] ~ 0.125 eV for w 0.02 eV. This value is smaller than the current 

cosmological bounds (see, e.g. [H]), but is within the sensitivity expected to be reached by 
combining data on weak lensing of galaxies by large scale structure with data from WMAP and 
PLANCK experiments (see, e.g. [12]). 

4.4 Specific Features of the Predictions for B{fi — >■ e + 7) 

Apart from discussing generic results for the AF and the AM model it is also interesting to have 
a closer look at particular points in the parameter space of the A4 models (with generic NLO 
corrections). In order to do so we use the analytic formula given in Section 3.2., eq. (|3.10p . for 
the branching ratio of the decay — )• e + 7 together with the results given in Table 1 and assume 
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that the coefficients of the 0(e) terms are real and have the same absolute value C > 0: 



\yl log + C e log log + .3 log ) . 

3 ymij \ J \mij V"^i / / 

(4.3) 

We do not fix the relative sign of these terms and allow for all eight combinations cri,2,3 = il- 
We choose for the mSUGRA parameters the values given in eq. (|4.ip . set e = 0.04, set yu = 1 
and take again best fit values for r and |Am^|. In Fig. [6l left panel, we plot the result for 
B(fi — )• e + 7) for C = 1.3 in the case of a NO light neutrino mass spectrum with respect to 

the effective Majorana mass niee, fnee — \2mi + m\ + AmQ|/3. We show only the two curves 
which correspond to the upper and lower bound that can be reached for the eight different 
combinations of i7i,2,3- As one can see, there exists the possibility of cancellations between the 
terms contributing to the branching ratio of the — t- e + 7 decay, so that the value of the 
latter can be strongly suppressed 1^. The value of at which the suppression takes place 
depends on the value of C. In Fig. [6l right panel, we show the corresponding plot for 10 
neutrino mass spectrum. We choose C = 1. In contrast to the case of NO spectrum, no strong 
suppression of — e + 7) is possible, because the term cx log(mi/m*) always dominates 
(see Fig. [21 right panel). This result holds for all values of the constant C from the interval 
0.1 < C < 6. Allowing for arbitrary relative phases between the different contributions in the 
right-hand side of eq. (|4.3p . we find that for a NO light neutrino mass spectrum the curve for 
ci,2,3 = +1 corresponds to an upper bound on B(/x — )• e + 7), whereas the curve for o"i^2,3 = — 1 
is an absolute lower bound with the exception of few points in the parameter space. For the 10 
spectrum, the bounds obtained for real coefficients are also upper and lower bounds in the case 
of arbitrary relative phases between the different terms in eq. (j4.3p . 

As mentioned earlier, the preceding analysis holds for an ^4 model with generic NLO correc- 
tions, as is the case of the AM model. In order to perform a similar analysis for the AF model 
V\ we replace ai log(mi/m*) with eai log(mi/m*) in eq. ()4.3p . We find that deep cancellations 
between the different contributions m. ^ e + ^) can occur in both the cases of NO and 10 
neutrino mass spectrum. For the 10 spectrum, the cancellations leading to a strong suppression 
of branching ratio — )• e -|- 7) take place for ruee around 0.09 eV for almost all values of C in 
the range considered, 0.1 < C < 6. 



B'{ii — )• e -I- 7) oc 



5 The fi — e Conversion and ii — )- 3£j Decay Rates 

We briefly discuss in this section the experimental constraints that can be imposed on the 
models from data on — e conversion and the decays ii — )• 31 j. In the mSUGRA scenario, these 
LFV processes are dominated by the contribution coming from the 7— penguin diagrams. As a 
consequence, for fi — e conversion, the following relation holds with a good approximation [22]: 

CT(,N ^ eN) . Ik^i^ = <^4^zyFm, ^ e + .) . (5-4) 

J- capt iZTT i capt 

In eq. (j5.4p Z is the proton number in the nucleus N, F{q^) is the nuclear form factor at 
momentum transfer q, -^e// is an effective atomic charge and Fcapt is the experimentally known 

^^Note that the value of _B(/i — >■ e + 7) will still be non-zero in general, because we expect corrections to the 
coefficients of the different logarithms of order . 

We remind the reader that in the AF model the coefficient of the logarithm log(mi/m,) is of order rather 
than e. 
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total muon capture rate. For 22Ti we have Zf.ff = 17.6, F{(^ = — m^) 0.54 and Fcapt = 
2.590 X lO^sec-i [33]. In the case of f|Al one finds Z^^f = 11.48, = -m^) 0.64 and 

Fcapt = 7.054 X 10^ sec~^ [H]. According to eq. ()5.4p . the /i — e conversion ratios in g^Ti and 
^gAl are given by: 

CR{fi ^Ti ^ e ^Ti) « 0.005B(/x ^ e + 7) , Cii(^ ?|A1 ^ e f|Al) « 0.0027S(// ^ e + 7) . 

(5.5) 

Future experimental searches for /x — e conversion can reach the sensitivity: CR{fi 22''^'^ ~^ 
e^Ti) PS 10-18 03], and ?|A1 ^ e gAl) ^ IQ-^^ [S]. The upper bound e + 7) < 

10"^^ which can be obtained in the MEG experiment would correspond to the following upper 
bounds on the // — e conversion ratios under discussion: CR{fi 22''^'^ ~^ ^ll^^i) < 5 x 10"^^ 
and CR{iJ,ilAl — ef3Al) < 2.7 x 10^^^. The latter could be probed by future experiments 
on yu — e conversion, which have higher prospective sensitivity. We see that for the mSUGRA 
points considered in Section 4, both the AF and AM models can be further constrained by 
the experiments on ^ — e conversion if the ;U — )• e + 7 decay will not be observed in the MEG 
experiment. 

In what concerns the decay of a charged lepton into three charged leptons, the branching 
ratio is approximately given by [22]: 

B{i^ ^ 3£,) « ^ l^log (^^^ - B{£, ^ i, + 7) . (5.6) 

The searches for r — )• ^ + 7, r — 3/i and r — )■ 3e decays at SuperB factories [35] will be sensitive 
to B{t — )• /X + 7),S(r — )• 3fi),B{T — )• 3e) > 10^^. Therefore, if in the experiments at SuperB 
factories it is found that -B(r — )• /i + 7) < 10~^, obtaining the upper limits B{t ^ 3//), B{t ^ 
3e) < 10"^ would not constrain further the models considered here. However, the observation 
of the T ^ 3fi decay with a branching ratio B{t ^ 3fi) > 10~^, combined with the upper limit 
B{t — )• /X + 7) < 10^^, or the observation of the r — )• /i + 7 decay with a branching ratio 
B(t ^ fi + ^) > 10^^, would rule out the models under discussion. 

The current limit on the ;U — )• 3e decay branching ratio is i3(;U — )• 3e) < 10^^^ ^45j. There 
are no plans at present to perform a new experimental search for the — ?• 3e decay with higher 
precision. 



6 Conclusions 

In this paper we have studied charged lepton flavour violating (LFV) radiative decays, fi — )• e+7, 
T ^ + j and T —7- e + 7 in a class of super symmetric models with three heavy RH Majorana 
neutrinos, in which the lepton (neutrino) mixing is predicted to leading order (LO) to be tri- 
bimaximal (TB). The light neutrino masses are generated via the type I see-saw mechanism. We 
work within the framework of the minimal supergravity (mSUGRA) scenario, which provides 
flavour universal boundary conditions at the scale of grand unification Mx ~ 2 x 10^^ GeV. In 
this scenario the flavour universal slepton masses, trilinear couplings and the gaugino masses 
are characterised by the parameters mo, ^0 and mi/2 respectively. The other free parameters 
are tan/3 and sign(|u). Flavour off-diagonal elements in the slepton mass matrices are generated 
through RG effects (from the scale Mx to the scale of heavy RH Majorana neutrino masses). 
The former induce the LFV decays /i— >-e + 7, r— 7';U + 7 and r — t- e + 7. In the present 
article the branching ratios B{ii — )• ij + 7) are calculated using the analytic approximations 
developed in p2H2^ . In this approach B{£i — )• £j + 7) depend only on the matrix of neutrino 
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Yukawa couplings defined in the flavour basis, on the three heavy RH Majorana neutrino 
masses, and on an "average" SUSY mass scale, ms (see eq. (|3.8p ). In the class of models 
considered, each of the three light Majorana neutrino masses rrii is directly related at LO to 
the corresponding heavy Majorana neutrino mass Mj, eq. (j2.15p . Both types of light neutrino 
mass spectrum (with normal ordering (NO) and inverted ordering (10)), are possible. The 
lightest neutrino masses, corresponding to the two types of spectrum, are constrained to lie in 
the intervals: mi ^ (3.8 ^ 7.0) x 10"^ eV (NO) and ms > 0.02 eV (10). The RG effects on the 
neutrino masses and mixing angles were not taken into account. As is well known, these effects 
are relatively small if the lightest neutrino mass satisfies min(mj) < 0.10 eV. For this reason we 
limited our analysis in the case of 10 spectrum to ms < 0.10 eV. 

The analytic estimates of the branching ratios B{ii — t- ij + 7) we have given were made 
for the case of generic next-to-leading order (NLO) corrections to the neutrino Yukawa matrix. 
The numerical results we presented, however, are obtained for two explicit realizations of the 
A4 models, those by Altarelli and Feruglio (AF) and by Altarelli and Meloni (AM) 0, 
respectively. In these models the ^4 symmetry is broken spontaneously at high energies and the 
symmetry breaking parameter e has a value in the interval 0.007 ^ e ^ 0.05. In the numerical 
calculations we set to e = 0.04. In this case the allowed range of tan /3 is 2 < tan f3 < 12. We 
used the value of tan/3 = 10 in the numerical calculations. 

The predictions for the £i — )■ £j ■ + 7 decay branching ratios, B{£i — >• £j + 7), in the AF and 
AM models are derived for one specific point in the mSUGRA parameter space lying in the 
stau co-annihilation region, which is compatible with direct bounds on sparticle masses and the 
requirement of explaining the amount of dark matter (DM) in the Universe: rriQ = 150 GeV, 
"^1/2 = 700 GeV and Aq = GeV. These values correspond to a bino-like LSP with a mass of 
approximately 280 GeV, which makes up the DM of the Universe, gluino masses of 1.6 TeV and 
squark masses between 1.1 TeV and 1.5 TeV. Having the indicated masses, the SUSY particles 
can be observed at LHC. Results for other points in the mSUGRA parameter space can be easily 
obtained by modifying the rescaling function present in the expressions for B{ii — )• ij + 7) (see 
eqs. dST]) and ([37711 ). 

We have found that in the case of NO light neutrino mass spectrum, both the AF and AM 
models predict B{fi — )• e + 7) > 10~^^ in practically all the parameter space considered (Figs. 
[3] and HI). The same conclusion is valid for the 10 mass spectrum in the case of the AM model, 
whereas for the AF model this result holds roughly in half of the parameter space of the model. 
Values of B{fi —7-6 + 7) ^ 10~^^ can be probed in the MEG experiment which is taking data at 
present. 

The predictions of the AF model for all the three branching ratios B{ii — )■ ij + 7) in the 
case of NO spectrum show a noticeable dependence on the value of the lightest neutrino mass. 
The dependence of i?(T —t- /i + 7) on min(mj) is particularly strong in the case of 10 spectrum. 
In contrast, B{fi — )• e + 7) and i?(T e + 7) in this case vary relatively little with min(mj). 
The predictions for B(ii — )• ij + 7) in the AM model do not exhibit significant dependence on 
min(mj). The branching ratios B{t — ?■ e + 7) and B{t —7-^ + 7) are always predicted to be 
below the sensitivity of the present and future planned experiments. We have shown, however, 
that if the SUSY particles are lighter, e.g. the LSP has a mass of 115 GeV, squarks are as 
light as 500 GeV and gluinos have masses 700 GeV, one can have B{n — )• e + 7) > 10~^^ 
and i?(r —)•// + 7) w 10~^ in the AF model with 10 spectrum (Fig. O left panel). A value of 
B{t —7- ;U + 7) ~ 10~^ requires the lightest neutrino mass to be 1713 ~ 0.02 eV. Sensitivity to 
such a value of B{t fj, + j) can be achieved, in principle, at a SuperB factory [H5] . 

We found that the dependence of the branching ratios B{ii — )• ij + 7) on the A4 flavour 
symmetry breaking parameter e in the case of NO spectrum differs from that in the case of 
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10 spectrum. For the NO spectrum, all the three branching ratios are predicted to scale as 
e'^. The normalized branching ratios B'{n — )• e + 7) ~ B(fj, — )• e + 7) and -B'(r — )• e + 7) = 
B{t — )• e + 7)/i?(r e + Vr + i^e), are found to be equal at LO, while B'{t fi + ^) = B{t — )• 
/i + 7)/i?(r — )• fi + i^r + ^^i) is typically larger by one order of magnitude. In the case of 10 
spectrum we have at LO B{n — )• e + 7) oc e^, B'{t — e + 7) oc e^. However, for e f« 0.04, this e 
suppression is compensated by relatively large values of the logarithm log (mi /m*). Therefore in 
a model with generic NLO corrections to the matrix of neutrino Yukawa couplings, — )• e + 7) 
and B'[t — )• e + 7) can have the same magnitude for the two types - NO and 10, of neutrino 
mass spectrum. The specific structure of the NLO contributions to the matrix of neutrino 
Yukawa couplings in the AF model leads to an additional suppression of the branching ratios 
B{^ — e + 7) and B'{t — )■ e + 7). As a consequence, B{^ — )• e + 7) and B'{t — )■ e + 7) in the 
case of 10 spectrum are predicted to be smaller in the AF model than in the AM model by up 
to two orders of magnitude. Since B'{t — )• ^ + 7) always receives a contribution of order e'^, 
it has the same magnitude for the NO and 10 spectrum, and in the case of 10 spectrum and 
for ms ~ 0.02 eV it is not affected by the indicated additional suppression. For the AF model 
thus i?'(r — )• /u + 7) can be by up to three orders of magnitude larger than — )• e + 7) and 
B'{t -)> e + 7). 

The type of dependence of the branching ratios on the symmetry breaking parameter e 
is in contrast to that found in the framework of an effective field theory approach in [15j, in 
which the setup of the AF model was used. This approach suggests that B{ii — )• ij + 7) scale 
generically as e^, and as in a specific super symmetric version of the theory. Moreover, the 
results found in the effective field theory approach do not show any dependence on the type of 
the light neutrino mass spectrum and on the value of the lightest neutrino mass. 

In [TU] we studied versions of the AF and AM models, in which the neutrino Yukawa couplings 
are suppressed, being proportional to e. As a consequence, the scale of the RH neutrino masses 
is lowered. Since the branching ratios of the charged LFV radiative decays studied in the present 
article are proportional to the Yukawa couplings to the forth power, suppressing the latter also 
efficiently suppresses the li £j + j decay branching ratios. Thus, for the mSUGRA points 
chosen in our numerical study and for most of the relevant parameter space of the version of A4 
models considered in pO] we expect i?(/U —)• e + 7) to be below the sensitivity of the on-going 
MEG experiment. 

We have shown that for specific values of the parameters in an ^4 model with generic NLO 
corrections, "accidental" cancellations between the different contributions in the — >• e+7 decay 
rate and a strong suppression of the branching ratio — )• e + 7) are possible in the case of 
NO neutrino mass spectrum. The same result holds in the AF model for both types - NO and 
10, of spectrum. No similar suppression is found to occur for the 10 spectrum in the ^4 model 
with generic NLO corrections. 

Our estimates of the predicted rate of /i — e conversion in the ^4 models considered show 
that future experiments can further constrain these models if the /U — )• e + 7 decay will not be 
observed in the MEG experiment. The observation at the SuperB factories of the r — t- 3/U decay 
with a branching ratio B{t — t- 3fi) > 10~^, combined with the upper limit B{t — t- ^ + 7) < 10^^, 
or the observation of the t ^ /jL + j decay with branching ratio -B(t — )• /Lt + 7) > 10~^, would rule 
out the models under discussion. If i?(r —)•// + 7) is found to satisfy -B(r —)•// + 7) < 10~^, 
the prospective sensitivity of SuperB factories to the decay modes r — )• 3// and r — )• 3e would 
not allow to obtain additional constraints on the parameter space of the A4 models from non- 
observation of the T ^ 3fj, and r — >■ 3e decays. 

The results of the MEG experiment which is taking data at present, and of the upcoming 
experiments at LHC can provide significant tests of, and can severely constrain, the class of 
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models predicting tri-bimaximal neutrino mixing, in general, and the Altarelli-Feruglio and 
Altarelli-Meloni models, in particular. 
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